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ABSTRACT 


A computer program is presented which solves the search 
planning problem for survivors at sea. The program is 
designed to utilize weather data as compiled by the United 
States Navy at its Fleet Numerical Weather Central, 


Monterey, California. 


An investigation is also made into the search criteria 
used by the United States Coast Guard in its planning 
procedures. Guidelines are given for the use of the Square 


search and the Sector search. 
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1. INTRODUCTION 


A. THE COAST GUARD AND SURVIVOR SEARCH 
Survivor search is defined as the process of planning 
and executing a search for survivors at sea, and the Coast 
Guard usually accomplishes this in the following steps: 
1. Initial Notification 
This stage involves gathering all available data 
about persons or craft in distress, including 
characteristics of the survivor craft, number of survivors, 
etc. 
~~ «(4 i Di Pos 
In this stage, the search planner puts together all the 
available data concerning the survivor and makes a judgement 
as to what the most likely position of the survivor was at 
the time the craft seased operating normally. 
ae) Ori Cal atio 
The survivor's suspected motion is plotted assuming 
drift in the open ocean, and a position is calculated for 
some future time, usually the planned beginning of search 
efforts. 
4, alculati f A b ar d 
An estimate is made of the accuracy of the 
survivor's position as calculated in step 3 and the accuracy 


of navigation of the search craft. An area, called the 


Search Area, is defined based on these errors. 
5. Allocation of Search Units 
A determination is made of the resources (ships, 
planes, boats, etc.) to be used in the search, and the way 
in which these resources will be used. 
6. Execution of the Search 
Searching is accomplished according to several 
established search patterns, most of which are different 
arrangements of a series of parallel! paths within the Search 
Area. 
7. Jermination of the Search 
Either the survivors are located or, after a certain 
amount of searching is completed without success, the search 
is terminated. This stage includes reports, and_ the 


evaluation of search effort. 


B. DOCTRINE 


The guidance for all these activities is contained in 
Coast Guard publication CG-308, the National Search and 
Rescue Manual, (Ref. 1.). Chapters 3 through 8 of this 
manual contain the instructions for planning and conducting 


a search for survivors at sea. 


CS RURPOSE 

This thesis accomplishes two things: 

1. Provides a computerized solution to step 3. of the 
search planning process, calculating survivor craft drift, 
and 

22 Investigates the criteria used for allocation of 


resources to a search, steps 4, 5, and 6. 


It is intended that the results be understood by the 


search planner, so that they may be of some use cio him. 


ll. THE SEARCH PLANNING PROGRAM 


A, OBJECTIVE 

The major objective in designing a computer program for 
search planning was to have ae fast, accurate means of 
obtaining an estimate of a survivor's oceanic position at 


the time of the beginning of a search. 


B. THE MANUAL METHOD 
In order to understand the program, it is firs 
necessary to understand how drift plotting is done manually. 
1. Example 
A fishing vessel is believed to have sunk at 
latitude 40-42.0N, and longitude 52-13.0W, at 09002 on 22 
June. The probable error of this position is estimated to 
be 30 miles. Although the standard measure of variability is 
the standard deviation, all error estimates made in_ the 
search planning process are in terms of probable error. It 
will be demonstrated that this error is defined by the 
statistical term, "probable error". Survivors are believed 
to be in rubber rafts. 


The proper Coast Guard authorities are notified at 


2100Z on 22 June, and a search plane has been dispatched. It 
is due to arrive on the scene at 0100Z 23 June. The 
aircraft's navigation error (again a probable error) is 
estimated to be 10 miles. 
2. Gathering Data 

The first task of the search planner is to estimate 
wind, wind current, and sea current at the position of the 
sinking. Wind is obtained from observations, weather maps, 
Navy Fleet weather services, or the U.S. Weather Bureau, 
whichever is available. Wind current is computed from the 
wind by means of a graph given in the National Search and 
Rescue Manual. Average Sea Current is obtained from ocean 
current charts or pilot charts. These sources for average 
sea current are compiled from observations made over several 
years. (Thus this data may not be accurate for the time and 
place of concern, but so far they have been the only data 
avallable to the search planner.) 

3. Orift Plotting 

a. Leeway 

Leeway is the motion of the survivor craft caused by 
wind. For the example case of a liferaft, two graphs are 
provided in the SAR Manual for computing leeway. One graph 
is for when a drogue is attached to the raft (a drogue is a 
type of sea anchor designed to retard liferaft leeway) and 
the other is for when a drogue is not attached. If it is not 


known whether a drogue is attached, both conditions are 


considered separately. 

b. Total Drift 

Once leeway, wind current, and average sea scurrent 
have been determined, they are added vectorially to produce 


total liferaft drift. 
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This calculation can be done for intervals of any 
length - one hour, three hours, twelve hours - but” the 
interval most commonly used is the interval between’ wind 
observations. That is if wind is available for 0900Z, 1200Z, 
1500Z, etc., the search planner would probably compute total 
drift in three hour’ intervals, from 0900Z on the 22nd to 
0100Z on the 23rd. Each time total drift is calculated 
(another probable error), 12.5% of that distance is retained 
as "drift error", and the drift errors are summed to obtain 
total drift error. 

Example: In this manner, the position of the example 
liferaft is calculated to be 40-58.0N, 51-55.0W, and the 
drift error is &.0 miles. 

c. Minimax Plotting 

lf it re unknown whether a drogue was’ in use, 
two plots are made, and the final ("datum") position is 


defined to be midway between the two positions calculated. 
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The drift error is a function of the two calculated drift 
errors and the distance between the two positions. This 
process of plotting two positions is called Minimax plotting 
in the SAR Manual. 
d. Error of Position 
Once the datum has been computed, jit is then 


necessary to find total probable error of position, c, where 


de = total drift error 
X = error in Initial positton of survivor craft 


Y = navigation error of search craft 
. 2 az ve 
c= de’ + X° +r Y 


For the example, c= \/16+100+900 =\/1016 = 31.9 miles. 


e. safety Factor 
To determine the size of the area to. be 
searched, a search "radius", R is found by multiplying c by 
a safety factor, depending on which search in a sequence is 
being planned. For searches one through five, the respective 
safety factors are 1.1, 1.6, 2.0, 2.3, and 2.5. The example 
of the sunken fishing boat is a first search, so R =1.1 x 
14.7 = 16.2 miles. 
The area to be searched, then, is a square with 


sides = 2xR, centered on datum: 


iss 


t 


ZR 


CO, 
DATUM na 


The standard measure of the effective search 


f. Sweep Width 


width of a search craft is called "sweep width". The value 
of the sweep width is obtained from a table in’ the Search 
and Rescue Manual, based on the type of target, search 
altitude, and meteorological visibility. The sweep width 
thus obtained is corrected for the effect of sea state by 
multiplying by a "“whitecap correction factor" which is 


tabulated for various wind speeds. 


C. THE COMPUTER PROGRAM 
1, Qbjective 
The computer program accomplishes the objective of 
providing an estimate of the survivor's position by 
automating the steps of the manual method of drift plotting. 
2. source of Weather Data 
!n choosing a source of weather data for the computer 
program, three criteria were considered: a. availability, b. 
computer access time, andc. accuracy. The final selection 
for data source was the Navy Fleet Numerical Weather 
Central, Monterey. 
a. Availability 


FNWC Weather data is compiled twice daily, at 


2 


0000Z and 1200Z, and iis directly available, either by data 
link or, for testing purposes, on magnetic tape. For use 
with this program, the data is stacked chronologically on a 
master tape. 

b. Computer Access Time 

With the exception of the mounting of the master 

weather tape, all weather data handling is’ eliminated. 
Access speed, then, is limited only by the speed of the tape 
transport. 

c. Accuracy 

Although the data used is intended for wide area 
coverage, its probable error is less than the probable error 
assumed when computing drift error (12.5%). Although there 
will be occasions when local weather reports are available 
and may be more accurate, the bulk of search planning is 
done under conditions of incomplete weather analysis’ and 
forecasts, and when gross estimates must be made. Since the 
weather analyses for a particular time period and 
geographical location, once assembled for this program, are 
fixed parameters, repeated computations of the same search 
problem will always produce identical answers. This permits 
a rigorous evaluation of the data, methods, and assumptions 
used in search planning by a systematic testing procedure. 
a. ~=—D Descri n 
The weather data are produced by FNWC for 3,969 


geographical points in the northern hemisphere. These data 
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points are intersections of a 63 xX 63 square grid 
superimposed on a polar’ stereographic map projection. Grid 
point (32,32) is at the north pole, and the vertical line 
(j=32) is coincidental with the 80 W and 100 E meridians. A 
sample projection is shown as Figure l, 

The polar stereographic projection provides a 
plotting reference which shows an entire hemisphere with the 
least possible distortion (Ref. 6), and the scale is true at 
60 degrees North Latitude. Conversion from latitude (L) and 
longitude (A) to grid coordinates is accomplished by the 


following formulas: 


N{- 


(ie Vest. 205 —— Cos (a +10 ) (1) 
(os Shee 


n\- 


i Li enh aa ane fee sin (x +10 ) (8 


/ +StnNL 





For the example, i = 38.008, and j = 19.0143. 


The data is stored in "fields" of 3,969 words 


(numbers), each word being some component of atmospheric or 


oceanic conditions at its corresponding (i,j) grid position. 


The fields required for direct Input to this program are: 
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Figure 1. Portion of Polar Stereographic Projection Showing 
FNWC (i,j) Coordinates 


i 





Oy CFR am U- component of surface current 


V Curr V- component of surface current 
U Marine U- component of surface wind 
V Marine V- component of surface wind 
where the U direction is parallel to the "i" or 


horizontal axis and the V direction is parallel to the "j" 
or vertical axis, in the FNWC grid system. 

The Current information as provided by FNWC is total 
surface current, including wind current. 

Fleet Numerical did not regularly produce U Marine 
and V Marine wind fields, so these fields are computed by a 
separate computer program provided by the Weather Facility 


using the following data fields as input: 


D 1000 Height of the 1000-mb pressure line 
T Air Air temperature at the surface 
T Sea Sea temperature at the surface 


D. CALCULATING DRIFT WITH THE COMPUTER 
1. \Leeway 
In addition to the liferaft drift tables given in 
the Search and Rescue manual, there are other ways’ to 
predict the leeway of a survivor craft. Most of these 


methods involve using a percent of the wind speed _ to 
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determine the leeway speed, and some deflection in degrees 
to determine the direction of leeway off the wind line. For 
example, a boat type may be known to drift in a direction 40 
degrees off the wind line, and at a speed equal to 4% of the 
wind. To accomodate these methods as well as the liferaft 
drift table from the SAR Manual, several options are 
available to the user of this computer program. 
2. Options for Leeway 

a. A function approximating the liferaft drift 

tables in the National SAR Manual has been provided; 


drifting speeds for with drogue and without drogue are both 





available. 

b. Leeway speed can also be specified as a percent 
of wind, in whole numbers, from 2% to 9%. 

c. For each of the leeway velocities in b., a 
deflection off the wind line of plus or minus some angle, in 
tens of degrees, from 0 to 90 degrees, may be specified. 

d. A minimum-maximum leeway speed option with the 
minimum being 1% of the wind, and the maximum being some 
whole number percent from 2% to 9% of wind speed, is 
available. 

e. No leeway at all (as in the case of a man in 
the water or some other almost completely submerged object) 
may be specified. 


3. Adjustment of Input Data 


a. Dates 
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All dates are converted from month and day to 
julian dates, to permit reference to weather data _ and 
adjustment of the program "clock". 

b. Latitude and Longitude 

All initial positions are converted to FNWC 
(i,j) coordinates. 

c. Leeway Options 

Leeway options are entered by means of a 
two-digit code, which Is converted to speed and direction of 
drift within the program. 

d. Weather 

Wind and current are provided in units of 
cm./sec. and are scaled in the program to nautical miles per 
hour. They are further scaled to decimal fractions of the 
FNWC grid distance. 

4, W h D 
a. Interpolation for Position 

For each position in (i,j) coordinates there are 
four neltghboring grid intersections for which weather’ is 
tabulated. Example: 


GRID POSITION: i=38.00787 


j=19.01413 
GRID POINTS: . (38,20) . (39,20) 
. (38,19) . (39,19) 


All weather data are interpolated linearly among 
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these four neighboring points. For example, if U Current at 
0000Z at grid points (i,j) is: 


(38,20) 50e0 (39,20) 


i 
_ 
(O 
© 
© 


40.0 


(38,19) 42.0 (39 bo 


then for {=38, j=19.01413, 


U Current = 42.0 + .01413 x (50-42) h 2012 

and for {=39, a similar computation results in U 
cCumrent = 40.12/17. 

U Current for the position is then, 

42.112 + .007387 x (=-1.985)> = 0275007. sec. 

b. Interpolation for Time 

For each hour in the computation, there are two 
observation times which are used to supply weather data: 
0000Z and 1200Z, one of which will be prior to the hour, and 
one which will be after the hour. A linear interpolation for 
position is done for each observation, and then these two 
values are interpolated linearly for time. For example, if 
the U Curr for 1200Z is calculated as above and found to be 
37.25 cm./sec., the interpolation for time is accomplished 
as follows: 

46.09 + 9/12 x (-8.84) = 39.46 cm./sec. 

c. The Clock 

A complete set of interpolations is done ffor 
each position and for each weather parameter at the start of 
the drift calculations, and each position is advanced for an 


hour's worth of drift. This interval is used as 
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compromise between accuracy requirements and computer time. 
If an interval much smaller than one hour were used, a large 
increase in computer time would be required. A_ problem 
which would take 20 seconds at one hour intervals would cost 
almost 20 minutes of computer time if one minute Intervals 
were used. 

After the one hour's worth of drift is computed, 
and the positions updated, the clock is advanced one hour, 
and the process is repeated over and over until the end of 
the problem. 

Ses ite P| 
Drift calculations are applied to the (i,j) 
coordinates of the survivor craft plot. That Is, the 
U-component of drift ts added to the "i'' coordinate of the 
survivor position and the V-component is added to the "j" 
coordinate. 
Example: The Interpolated values of current = and 


wind for 0900Z on 22 June are 


Parameter cm./sec. miles/hour 
U Curr + 39.460 + .766 
V Curr + 36.128 + .702 
U Marine +250.78 + &,87 
V Martine -850.25 -16.55 


the wind speed then equals \/4.61'+ 16.55" #17.25 
nautical miles/hour. For a liferaft without drogue, this 


wind speed produces a leeway of 0.94 miles per hour as 
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follows: 

Leeway = .226 x\f 17.25 = 0.94 

This formula is an approximation to the liferaft 
drift table given in the Search and Rescue Manual. 

Stnce arctan (4%.87/-16.55) equals 163.6 degrees, 
this angle is a measure of the wind angle (measured 


clockwise from the vertical (Cin the FNWC system): 


If leeway deflection is specified by the user, it is 
applied to this wind angle, before the U and V components of 
leeway are calculated. The reduction of Leeway and wind 
angle to U and V components is as follows: 


0.94 x sin(163.6°) 


0.266 


U Leeway 


V Leeway 0.94 x cos(163.6°) -0.902 


The leeway and current are applied to the position 
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as follows: 


U V 
Leeway . 266 ie JZ 
Current _+/66 ested ads 
Total 1052 -.200 knots 


These values are then divided by the distance 


between gridpoints, at this latitude = 235.37 miles: 


00439 =, 00035 
+ position 38,00787 19,01413 
New Position 38.01226 19.01328 for 1000Z 22 June. 


This process is repeated hour by hour until the last 
datum has been calculated, drift plotting ceases, and 
calculations for position error and Search Radius begin. 

6. Position Error 
Just as in the manual method, the _ initial 
position error of the survivor craft is computed as follows: 

x = distance of initial (or distress) position from 
last known position. 

A = error of last known position 

B = navigation error of survivor craft as a percent 
of distance traveled 

The initial position error , X, is 

X = A+B x 

Since Y is the navigation error of the searcher, 


and de is the drift error, the Total Probable Error, c, of 


“7 


Position Is 

c = \/X" + Y+ de? 

des a Fact 

c, Total Probable Error of Position Is multiplied by 
the safety factor to find the Search Radius as prescribed in 
the Search and Rescue Manual. 

8. Sweep Width 
A numerical approximation to the data In 
Richardson's Sweep Width tables (Ref. &) may be seen to be: 
W = -3.34377 + 1.25887 log (1) + 2.55442 log (V) + 
~01023 h 
where 1 is boat length 
V ts vistbility in miles, and 
h ts search altitude In thousands of feet. 

[In addition, corrections’ to Sweep Width for 
whitecaps based on wind speed (v) and cloud cover (C) can be 
approximated by: 

Fl = .886 + .298X - .116 X* (whitecaps, X = v/10) 
F2 = 1.131 - .472C + .152C* (cloud cover) 

These correction factors are multiplied by W_ to 
obtain final Sweep Width. This Sweep Width function Is 
applicable to small boats only. No Sweep Width for large 


ships or for liferafts Is available In the program. 


E. VERSIONS OF THE PROGRAM 


The Search Planning Program is written in FORTRAN IV for 
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use ona remote computer terminal. A second version is 
written for use onthe operations computer at Fleet 
Numerical Weather Facility. No remote terminal operation is 
available there, so these two versions are different In 
several ways: 

1, jneut-Output 

a. The time sharing version Its written for entry of 
data by the search planner onto a remote console, and the 
results are printed out immediately on the same_ console. 
This permits the planner to control the program during 
execution, and to make intermediate decisions concerning the 
problem. An example of the output of the time-sharing 
version is shown on the next two pages. Note that items 
typed in by the user are shown in lower case, and the 
replies typed on the console by the computer are in upper 
case. The final output Includes sweep width, total probable 
error of position, search radius, and calculated wind. 

b. The batch-processing version at FNWC differs 
from this In that Itt requires cards to be punched for input, 
and the output Is _ printed In message form, ready for 
transmission back to the requestor. 

2. Features Omitted 

For the operational version of the program at Fleet 
Numerical Weather Facility the following Items of output 
shown on the’ sample time-sharing printout were not made 


available: 
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$ sar 
EXECUTIOf SEGINS... 
MONTERFY SFARCH PLANNING PROGRAM 
ENTER: CEEXAY COU aie OA ew | EME ot RiP fice 


libohiteratd 
HOW MANY STARTING POINTS WILL PBF FNTFEREN? 
1 


LIST CHROMOLOGTCALIY: STARTING TIMES AND POS hi VO"\ Se ae 
O70930Z MAR, 3822. SNAGae Gnas 
220900zZ jun,4N42.0n,05215. Ow 
40 42.0 od 130 CHECK ah \) 
CARO Kael ero OR on 


y 

ENTER NUMBER OFVTIMES DATUM [Sere er CO PUT a tie 
2 

BNTER THESE DIGS, EXs@20900 MAR? 035 1,0 50gMAP 

ZOOINO stn, 250700 1un 

SEARCH 1 7SRRACK=e eS DATUM: 40 50.21 St Loves Gib 
SEARCH 2 RACK 661, Aue: LO 35.20 51 44.04 
ENTER EAST RNOWN POSTUION: EX: 3522.0N,07653.00m 
4100.0n,05400. 0w 

Ss ERA ||) 

Pte mOAl, “Ceesnie, OR (S74) RCRAPT? CENTERS LOR» 
1 

PV ER SEARCH UNTY FIX ERROR, Y., EXeoo0 

>. 0 

IDI CARED SOISTRESS CRAEL POS ITLIGM ERROR: X: 

a 15 OOM r+ OCIS 2) TIMES Hist. FROM BAST RO OS ie. 
GOR RmChoOnoe C1 ¥PE a0, 172 (9 OR ai?) 

1 

ENTER PARAMETER(1): 

2020 

[Se ACH WATUM TIME AYO TEFERE MT SPARC (ire 7 55 anne: 


no 

eee NUMER CORRES POND is TO.) Teh 
** J] 0 SEARCH: 

ir 

wear) SEARCH. 
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0 
THF FOLLOWING ARF THE DATUM POINTS FOR THE 
ls OO SFEARGH: 
bN-50.2N 5$1-52.0W POS. ERROR: 4128, RAPIUS: 6000 i ieee 
WIND 34.4KTS FROM 267, 
SWEFP WIDTH? Y OR N 


SY 
FNTER SEARCH Atm ee Cex. 2800...) 
S 10602 
ENTE Ray 1 SB Ii NM hiss 2CF 2100) 
> 20Re 
ENTER BOAT LENGTH IN FFET (EX:36.0) 
Sw. 
FNTER CLOUD COVER !N FRACTION (EX: 0.9) 
Or 
SWEEP WIDTH= LO MIbLES.. 


THE FOLLOWING ARE THE DATUM POINTS FOR THE 
oe OER CH: 
4O0-53.3N 51-44 .0W POS. ERRORS Wage aul Se 0.0 MipbESss 
WIND 34.,4KTS FROM 267, 
GOON LUCK ON YOUR SEARCH! LET ME KNOW TF YOU FIND 
THE Lt Pewee 


1HC0021 STOP 0 
Rs T=1.00/2.87 13.11.18 


Note: Inputs marked with ‘>! 
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a. Printout of intermediate track positions (SEARCH 
1, TRACK 1, etc.) is not produced. 
b. On scene Wind 
Since the wind printout as shown in’ the sample 
time-sharing output Its only the wind velocity and direction 
at the beginning of the search, and may not be 
representative of the average weather conditions throughout 
the search, this value is not included. 
c. Sweep Width 
Since the the computation of sweep width 
requires data which may not be immediately available to the 
search planner at the time of his request, and since the 
formula applies only to small boats, this feature was not 
included in the Fleet Numerical version. 
d, Leeway 
For the Fleet Numerical version, leeway’ is 
calculated as described in D.2. The time sharing version 
uses a square root function for all leeway. If N is the 
percentage drift specified (1% to 9%), the leeway function 
in the time sharing version is 
LeewAyY = “Ee an (3) 
This function produces leeway which is equal to 
NZ of the wind for wind speeds of 16.0 knots. For other 
wind velocities the function follows the general contour of 
the liferaft leeway (without drogue) table in the Search and 


Rescue Manual. 
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Pll. SEARCH THEORY 


In part II! the method of estimating the position of the 
search object and the rules for determining the size of the 
search area were discussed. This is only part of the search 
planning process; this section will consider how resources 
are assigned to the search, and how the effectiveness of the 
search is measured, 

In a search pattern employing parallel sweeps, the 
distance between these sweeps is called track spacing, S. 
The Search and Rescue Manual recommends’ that for urgent 
cases the track spacing used be equal to the Sweep Width, W, 
and for less urgent cases it be equal to twice the Sweep 
Width. Once the value of track spacing has been decided 
upon, the amount of searching required (in miles’ flown) L 


can be computed from S and the area A to be searched: 


L = A/S 


The measure of effectiveness of a search which is given 


in the Search and Rescue Manual ts called Probability of 


Detection, and it is presented as a series of curves plotted 


against "coverage factor" or W/S. These curves are called 
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"First search", "second search", etc., probabilities. If p 
is the probability of detection in the first search curve, 


the probabilities 2,..n..5 are given by 


p =1-(1- p) 


These curves do not take into account the size of the 
area searched, however, so that a search can be conducted In 
avery small area using a high coverage factor which will 
yield a very high probability of detection according to the 
graph. Common sense Indicates that to confine the search to 
a very small area will probably not result in the detection 
of the target. Nowhere does the Search and Rescue Manual 
define the probability that the target is in the search 
area, 

The absence of the treatment of this seemingly critical 
factor led the author to a review of the field of search 
theory as it applies to survivor search. 

The basis for most of the search planning criteria in 
the National Search and Rescue Manual was found in Koopman's 
(1944) work, Search and Screening (Ref. 2). This treatise 
developed the concept of sweep width, and found the 


probability of detection for various situations. 


A. DETECTION LAWS 


5 ig l n ous Prob Li of De tion 
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If Yis defined so that the probability of detecting 
a target ina short period of time dt Is Ydt, then the 
quantity Y Is called the instantaneous probability density 
of detection. When the searching is done during a time t 
under constant conditions, the probability p(t) of detection 


is given by 
~_yC 
(t) = —_< 
i (4) 


Search and Screening developed this quantity for 
various situations. 
a. Definite Range Law 
This law assumes that there is a certain range R 
from the observer beyond which it is impossible to detect 
the target; inside this range detection is certain. 


Instantaneous probability of detection by this law is 


Yr) = oo r<R 
(5) 
Y(r) = O r>R 
b. Inverse Cube Law 


The inverse cube law is the result of making the 
following assumptions: 

(1) The observer Is at an altitude h above the 
ocean, 

(2) The observer detects the target by seeing 
its wake. 


(3) The Instantaneous probability of detection 
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is proportional to the solid angle subtended at the point of 
observation by the wake. 
The calculation of the solid angle is shown in 


figure 2. 





Figure 2. Calculation of Solid Angle 


Solid Angle ex L . 
x= t/s ) IZ aie 
pur c/a = h/s 


thus ¢ =ah/s 4 
2 
Solid Angle = albh/s = abh /(r +n ) 


Since Y was assumed to be proportional to the 

solid angle, ea 
~ (heere) 

where k included all factors such as 


meteorological visibility, etc. 


31 


Assuming further that in the majority of cases r 


is much larger than h, then 
ee cll 
| ees » 3 





Dee ateral Range 

When the observer and the target are on their 
straight courses at constant speeds for a_ long time before 
and after their closest point of approach the probability of 
detection may be shown to be a function only of the distance 
at closest point of approach. This distance is called the 
lateral range, and is denoted as x. This probability of 
detection, p(x), plotted against lateral range, is called 
the lateral range curve. 

Twice the area under the lateral range curve Is 


defined as effective search width (or Sweep Width), 
too 


W= 4 andes (6) 
O 


Koopman shows that: 


"The effective search width is twice’ the 
range of a definite range law of detection 
which ts equivalent to the given law of 
detection in the sense that each of the two 
laws detects the same number of uniformly 
distributed targets of identical velocity." 


Thus, for the definite range law of detection p(x) =# 
l1 if xg R, and p(x) = 0 if x > R, so that W = 2R. To find W 


for the inverse cube law, first p(x) may be shown to be 
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2kh 
yy 
p(x) = i —- ¢ (7) 


—— 


where w is searcher's speed. 





From (CG), oan 
-f 30 = to xe 
Yn yma Sk 
_ eo (8) 
= > | 2tkh oe 
GJ 
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To find the probability of detecting a target known 
to be in an area A, where both searcher and target are 
moving over the ocean in complicated paths, the situation 
called random search may be considered. This requires’ the 
following assumptions: 

a. The target's position is uniformly distributed 
in A, given it has not been detected 

b. The observer's path is random in A in the sense 
that it can be thought of as having its different (not toc 
near) portions placed independently of one another in A. 

c. W is small with respect to the area A. 

If the searcher's path, of length L, can be 
considered to be composed of n separate and equal portions 
of length L/n, then the probability that the target be 


detected in any one path is WL/nA. The probability that the 


target not be detected at all, then is 


WL . 
p= | el 
(10) 
and for large n, 
~WL/A 
p= f-e (11) 


This is the formula of random search. 
4. Parallel Sweeps 
The most common method of non-random searching is 

one which uses a_ series of straight parallel sweeps a 
distance S (or track spacing) apart, where the position of 
the target is fixed with relation to the sweeps. Let the 
target be uniformly distributed in an area of size A. 
Search and Screening gives the following result for P(L), 
the probability of detecting with path length L=A/S. 

(1) Definite Range Law 

Since detection will surely occur if the target 
is within the detection range W/2, then’ the probability of 


detection is equal to the coverage factor W/S. 


Thus, 
p = WL /A i_< A/w 
(12) 
a { 1 > A/w 
(2) Inverse Cube Law 
T\yam alin ‘oa 
p= ert (“")- ert (5 A (13) 


(3) Uniform Random Search 


From (141), a ee 
p= { — eC 
h S=L/A, | 
or when / Le 
ya=tl —- <€ (14) 


The three probabilities using parallel sweeps 


are shown in Figure 3. 


DEFINITE RANGE ce 


{1 |j-—-—--—-—— — Ae, 
ee SEARCH 
Pi(s) WA 
: INVERSE CUBE LAW 


W /S 


Figure 3. Probability of Detection 


B. DISTRIBUTION OF SEARCH EFFORT 
l. General 
Another result in Search and Screening involves the 


definition of search effort density as a_ function of (x,y) 
in A, 


if B is the area of a subregion in A, and L is the 
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length of the observer's path in B, then the expression L /B 
is defined to be the observer's track density in B. As B 
becomes very small (approaches a point (x,y)) the search 
density @/ (x,y) is the limit of WL /B, or W1(x,y), where 


[tos yraedy iy ae (15) 
d 


an 
ee (16) 
Jf oe yJdady = Wh” & 
4 
Thus if “(x,y) is applied at (x,y), the probability 


of detecting a target located there is 
¥ ~ €(x%y) 
p — 
(17) 
Given the search density function (x,y) and the 
target's distribution p(x,y), the probability of detecting 
the target is | = PiOay) 
ple) = ff(e%y)Ci— & )dady (18) 
A 
The Y (x,y) which maximizes P( Y ) is shown by 


Koopman to be 


_ P(K)Y) 

P(X j)y) = log be 
(19) 

when (x,y) is in A and 
G(x,y) = O (20) 


/ 
when (x,y) is not InA , 


where AR denotes that portion of A in which” the 
expression (/49) is non-negative. The value of b is 


determined by using equation (16). 
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2. Normal ‘g Di ibution 
if the target's distribution is a bivariate circular 
cartesian 


distribution about a point(0,0) in a 








norma! 
coordinate system, a 
pity) = oh, S i ; Ga (aay 
(21) 
then the search effort density which maximizes P(@) 
is 
2 L 
OC —— aa< cs 
(22) 
and 
QD (x,y) cial) Gee 
(23) 
Thus letting a be the radius fr at which p(x,y) 
assumes the critical value b, " 
Aw=To ) b= ge.e 2 (24) 
where 
Se to" E] 4 C25), 
The form of this equation is an inverted parabola of 


with a maximum value of 


revolution with its center at (0,0) 


a ihe 
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IV. THE UNIFORM COVERAGE SQUARE SEARCH 


Although ( 19 ) gives the "optimum" search effort 
distribution, this optimum distribution is difficult to 
achieve in practice. Flying in a large circular area with a 
continuously varying effort is no simple navigational task. 
Much easier is the common practice of (1) defining a square 
area with its center at the center of the target's 
distribution and (2) conducting a regular search of uniform 
density throughout this area. If the target is not found, 
then (3) the area is expanded and the search repeated. 
Although a single search of this type does not conform to 
the parabolic optimum search effort distribution, a series 
of them, each larger than the one before, may be a very good 
approximation. 

The most Important advantage of this plan is its ease 
of navigation, allowing the pilot(s) free time to look for 


the target. 


A. THE PROBLEM 
The problem was to determine the optimum size of a 
square area within which a unltform search with parallel! 


sweeps is to be conducted. 
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B. TERMS 
The terms and definitions in the Search and Rescue 
Manual differ from those in Search and Screening. In the 
Search and Rescue Manual, probability of detection is used 
only for the probability of detection with parallel sweeps, 
(given a uniform target distribution within the search 
area). In order to avoid confusion, the following terms are 
defined: 
P - Probability of Success - The probability that on a 
given search the target will 


be found. 


p - Location Probability 2 The probability that a target 
is within a specified area. 


d - Probability of Detection - The probability of 
detecting 3 uniformly 
distributed target given 


parallel sweeps. 


W=- Sweep Width = The effective search width of the 
searcher; the area under” the 
lateral range curve. 


S - Track Spacing - The mean distance between adjacent 
parallel sweeps of the 
searcher. 


L - Path Length - The total number of miles flown by the 
searcher within the search 
area. 


c - Probable Error of Position - Probable error of the 
target's position from datum; 
given a normal distribution, 
c=1.18 x standard deviation of 
the target's distribution. 


R - Search Radius - The search radius tis defined = as 
one-half the side of a square 


bee) 


search area. 


C. THE SOLUTION 
1. The Objective 
The objective is to maximize the Probability of 
Success P, where the detection function for a uniformly 
distributed target known to be in the area ts’ defined as 


d(W,L,R,x,y). P is expressed as 
+- oe 






P(W,L,¢,R) = pias y lind ( Weill Sey ) dxdy 


For a untform search, the probability of detection is 
constant and equal to d(W,L,R) Inside the search area A, and 
is zero outside, (given unitform visibility and other 


factors), and so d(W,L,R) can come out of the integral: 


P(w,t,c,R) = d(w,L,R) p (Gy) dady 


It is convenient to choose x and y as distances ina 
cartesian coordinate system with the target's most likely 
position (datum) at (0,0). 

2. Calculation of Individual Terms 


Based on the inverse cube law for parallel sweeps, 
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Koopman derives the following expression for probability of 
detection, from (13): 
h ( zz WwW 
jie Cir je 
Z 
where S = A/L = (2R )/L. 
This is equivalent to the "First Search" probability of 
detection given in the National Search and Rescue Manual. 
The location probability, p, is given by the formula 


' for the circular bivariate normal distribution: 


x, ark 
Bm 
_ Ya ha 
p (RI Me € Ax dy 
Kolin ~R 


Where R defines a square centered at (0,0) with sides 
=2R. 

Figure & is a graph of this function, along with a 
similar curve showing the cumulative probability for a 
circle with radius R. Since the square area is actually 
larger (it is circumscribed about the circle) its 
probability for any radius R' is greater. However, if these 
cumulative probabilities are plotted against the areg of 
the search, the probability for a square will be slightly 
less than for a circle (see figure 5). Although there is 
some loss of efficiency in searching a square instead of a 


circle, this loss is very minimal, about 1%. 
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NOTE: Although the standard deviation ts the commonly 
accepted measure of variability, the National Search and 
Rescue Manual uses the term, "Total Probable Error of 
Position". Its value is computed ustng estimates for drift 
error, target's position error, and search craft navigation 
error. These errors are squared and added, and then a 
square root is taken to produce Total Probable Error. This 
term i[s not specifically defined in that publication, so 
that one of the objectives of this solution was to 
demonstrate that it is equivalent to the statistically 
well-defined term ‘Probable Error’, which Is a median error: 
that is, half of all errors have absolute values greater 
than that error, and half have absolute values which are 
smaller. The absence of this definition is a gap tn the 
present search planning process, because planners are called 
upon to make estimates of errors without knowing how these 
errors are defined. 
3. Ihe General solution 

To determine the optimum size for a square search, 
given a certain amount of searching miles available (L), 
sweep width (W), and probability distribution (p(x,y)), a 
straightforward maximization method its employed in_ the 
following sequence: 

a. Determine total effort, E. 


E=WxL 


Kh 


b. Choose some arbitrary initial search radius, x. 
c. Determine the search area, A, for this x: 
A = 4 x * 


d. Determine effort density for a uniform search: 


WL € Y.) 
A s 
e, Determine probability of detection d from (/3). 


f. Determine the probability p from (2!) with R = 


g. Determine the probability of success, 
P(x) = p(x)d(x). 
h. Increase x by Ax and repeat calculations. If 


P(x) increases, continue adding Ax until a loca? maximum Its 


encountered. If P(x) decreases initially, reduce x by 2Ax 
and continue subtracting Ax until a local maximum is 
reached. 


Since a cumulative probability distribution is 
continuous and has a slope throughout its range of > 0, and 
since the probability of detection function d(x) is 
continuous with a slope of < 0, the function P(x) = p(x)d(x) 
can be shown to be unimodal, and so a local maximum in P(x) 
will be an absolute maximum. 

4. Applica n 

The first use of this method was to heck the 

National Search and Rescue Manual's recommended search 


planning guidelines to see if they produce square searches 


of optimum sizes. These guidelines are given in the form of 
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'Safety Factors’, or search radii, for repeated searches, as 


multiples of Total Probable Error of Position, e.g., 


Search Safety Factor 
1 1 
2 1.6 
3 2.0 
ly Deo 
5 ee, 


As far as effort density is concerned, there are two 
recommendations given in the SAR Manual. One is for a track 
spacing equal to sweep width, for urgent cases, and the 
other for a track spacing equal to twice the sweep width, 
for less urgent searches. Both of these recommendations 
were tested. | 

The approach used was: 

a. Set some arbitrary value for W 

b. Set S=W so that WL/A. =1. 

Gr Find R from "safety factor'' table for the 
appropriate search (starting at 1). 


d. Find path length, L: 


4R 


ieee Ss 


e. Beginning at x = R, and with W and L above, use 


the general solution process described in paragraph 3. to 
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find x', the optimum search radius for that path length. 

f. Assuming that the optimal search Is executed 
unsuccessfully, reevaluate the target's distribution p(x) 
(see Para. 5.) and repeat the process from c. above. After 
five successive searches have been computed, stop_ the 
computation. 

zg. Compare the resulting five optimum sizes with 
the sizes recommended in the Safety Factor table from the 
Search and Rescue Manual. 

5. Bayesian Analysis 

Bayes’ Law of Probability says that an estimate of 
the probability of a hypothesis can and should be revised on 
the basis of the outcome of a_ related event. In the case of 
search theory, the hypothesis is that the target is inside a 
given area. A search its conducted, and the event is either 
the detection or the non-detection of the target. This 
event should revise the initial estimate of the probability 
that the target was in the area. Of course, if the event is 
detection, the new probability that the target is in the 
area is 1. If the event is non-detection, then Bayes' 
formula must be used to calculate the new probability. 

1. [£xample 

The probability that a target is in a certain search 
area has been determined to be 50%, and a search is 
conducted in this area with a unii@ara probability of 


detection of 50%. The resulting probability of success is .5 
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For the use of Bayes' formula, the following terms 


are defined: 


A The hypothesis that the target is in the 
area. 

A The hypothesis that the target is not in 
the area. 

D The event of detection. 

D The event of non-detection. 

Then 
| Ave (pla) 
old ld) = ) E 


pA). PLUIA) pla) plpdia) 


When the target Is detected: 


plald) = a eal’. 
O.5 x Ons. ¢O.5 «O10 
as expected. lf the target is not detected, the new 


probability would be 


o(41B) = p(A) + p(DIA) 
p(A)- e(D1A) + e(A)- P(DIA) 
S OCS a eal | 
(Lo eae ~~ ne 
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Whereas the target had a 50% a priori probability of 
being in the search area, its new, or ae posteriori 
probability of being in the area is 33.3%. 

For the solution in paragraph 4&., the entire 
probability distribution of the search object was 
reevaluated between successive searches. This was done for 
p(x) inside the search area by 


/ = p(x) i-d) 
ae (1 ~ p(R>d) Ae 


and for p(x) outside the search area by 


PG Vict la 
(1 — ptR)d ) (27) 


where R is the radius of the search conducted. 


In order to illustrate the general form of the 
probability distributions between searches, a function was 
used which is similar to a density function. This function, 


f(x) was defined to be 
d (p(x)) 
F(x) = ( aX ee) 


where x is the measure of radius, or again, one-half 


(28) 


of one side of a square, and 8xdx is the differential of 
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area (Figure 6). 
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Figure 6. Derivation of Density Function for Square 
Area 

Curve l1 in Figure 7 is the density function before 
any searching is done, as derived by (28) from (21). Curves 
2 through 5 are these functions after successive searches. 

Figure 8 is asimilar set of curves showing the 
cumulative probability distributions. 

6. Results 
This method (from &.) gives the following optimum 


searches: 


Search Optimum Radius Safety Factor 
1 1.08 ere 
2 1.64 1.6 
3 7.07 200 


50 


19 


~ 
re 


3 


~~ 
J - 


O96 


Density 


03% 


Probability 


Figure 7, 


cy 
= 
ut 
Lee 
i] 
cS 
C3 
KI 
OT 


Jus O1c 


R/f/c 


Probability Density with Application of Bayes' 
Between Successive Optimum Square Searches. 
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Since these results were very close to the ‘Safety 
Factors' listed above, it was concluded that (1) the safety 
factors were originally good calculations, and (2) they were 
calculated using the mathematical definition of probable 
error. 

It should be noted, however, that these safety 
factors are close to optimum search radii only for’ the 
condition where track spacing equals sweep. width. The 
second recommendation, track spacing equals twice sweep 
width, does not produce anywhere near optimum search areas. 
When a track spacing S' equal to twice the sweep width is 
used, the new path length L' is given by 

4R° 4er* | 


= = =o 
ies. 2W 7s 








For this new path length, optimum search sizes were 


found to be: 


Search Optimum Radius 
1 D8 5 
Z F20 
3 166 
P le 35 
> 2.20 
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When there Is less search effort avallable the sjze of 
the search area as well as the effort density should be 
reduced. The SAR Manual gives no guidance as to how much the 
area should be reduced with a given amount of available 
effort. Conversely, there is no indication as to how excess 
search resources should be used, if such a condition should 
exist. 

What Is needed is some way for the search planner to 
determine the optimum search configuration directly from the 


amount of search effort available. 


D. THE GENERAL CASE 

Assume that a series of optimum square searches, each 
larger than the preceding one, Its a good approximation of an 
optimum search effort distribution. Consider all the 
searches of the series combined Into one search effort, so 
that the search effort density @(r) assumes a step function 


which Is close to the ftdeal effort density. 


AcTUAL SEARCHES 


optimum EFFORT Dist. 






The radius r of the combIned search Is equal to the 


54 


radius of the last (and the largest) individual search in 
the sequence. 

if the search had been a circular one, and if the effort 
density had conformed to the optimum, the expression for 


this radius would be 


4cq*r 


= 


i 
4 





co 


But since the search is composed of a small number of 
Square searches, this same formula is not expected to apply. 
There should, however, be some way to approximate a function 
which will define the optimum radius based on available 
effort and the target's distribution. The simplest way to 
do this is to compute total search effort for a few optimum 
square searches, plot the radius of each of the searches 
against the cumulative search effort, and attempt to smooth 
a curve through them. 

This was done for the ten searches (five at W = S and 
five at W = 2S) calculated above, and the resulting curve is 
shown in Figure 9. Just as radius is in terms’ of probable 
error, so are the factors W and L of total effort 
"normalized' by dividing each by probable error. Thus 
"normalized' effort is a multiple of the square of probable 
error. 


it should be cautioned that the curve is only. an 
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Figure 9. 


Optimum Size of Square Search based on 
Cumulative Normalized Effort Factor. 
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approximation, and will not be valid if the effort available 
for a single search is very large. Thus if a search is 
computed which shows a coverage factor (W/S) greater than 
1.8, then the results of the curve will not give an optimum 
search, and some other method will be required. Again, the 
effort plotted is cumulative, so that all the effort exerted 
in earlier searches must be added to the effort available 


for the present search before entering the table. 


xample: One 150-knot aircraft is available for six 
hours of searching; the total probable error of position has 
been calculated to be 31.9 miles, and W is estimated at 10 
miles. 


Total effort then is found as follows: 
E= 150kts X 6hrs X 1lOmiles(W) 
= 9000 square miles 


In terms of probable error, (c=31.9) 


E(normalized) = 9000/(31.9 X 31.9) = 8.9 


From the graph, r = 1.4, 


1.4 X c = 44.6 miles, or 89.2 miles on a side, and 


* 
" 


4X RXR /L = 8.75 miles. 


2) 
i) 
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lf the target is not found on the first search, and a 
second is to be conducted, the amount of available effort is 
calculated for the second search. The normalized effort for 
this search must be combined with all previously expended 
effort before entering the graph. Thus if E(normalized) for 
the second search is calculated to be 10.0, the graph should 
be entered with a figure of 10.0+8.9 or 18.9 to find the 

radius of the second’ search. Each search area will be 
larger than the preceding one. 

This decision curve should prove to be amore flexible 
tool for search planning. than the Safety Factor table, since 
it provides a continuous range of search radii for different 


amounts of effort available. 
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V. THE SECTOR SEARCH 


There is one search pattern in general use which does 
effect a greater concentration of effort at the center of 
the search area, and that is the sector search. This search 
pattern is executed by a series of sweeps crossing over the 
center of the search area at various angles. These sweeps 
are connected at their endpoints by chords flown near the 
perimeter of the area. This pattern covers a circular area 
and, perfectly navigated, produces a track as shown in 


Figure 10. 


A. SEARCH EFFORT DISTRIBUTION 
The total effort expended in a sector search can be 


expressed by the integral 


E=f o(r) 2urdr 


© 
where R is search radius 


r is the distance from the center of the 


search 
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Figure 10. The Sector Search 
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@(r) is the search effort density at some 

distance r. 
This density function @(r) is simply W times’ the 
number of paths which intersect a circle of radius r, 
multiplied by dr to give total path length between_r and 


r+dr, and then divided by the area between r and r+dr. Thus 





Wn dar evi rs 
gir) = l > Fae (19) 
2Tr er yatta 
Then, oe 
E=WL = Wn PM ae = Wn (30) 
mo 
8) 
and, 
L=nR 


as expected. 
Here n is the number of ‘'spokes' produced, or twice the 
number of complete sweeps across the circle. This density 


function is illustrated in Figure 1l. 


B. PATTERN ERROR 

There will always be some error in flying any search 
pattern, which should be taken into account for the sector 
search. The navigation error of the searcher’ has already 
been taken into account when computing the Total Probable 


Error of Position. The pattern error, however, can be 
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Figure ll. Search Effort Density in a Perfectly 
Navigated Sector Search. 
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different from this navigation error. Whereas navigation 

error considers only the searcher's ability to find the 
datum point in order to define the search grea, pattern 
error, as it is we here, describes his flying error within 
the area. If, once the searcher finds what he thinks to be 
datum, he drops a smoke float or some other marking device, 
then his eaerern error will be much smaller than his 
navigation error. Conversely, if when arriving in the 
search area he stops’ refering to his navigation aids and 
relies thenceforth on dead reckoning only, his pattern error 
will be greater than his navigation error. 

Although by its nature, a flying error is’ bivariate 
(like the error of target position) a simplification was 
made. Since the pilot is flying long sweeps’ which have 
their endpoints relatively far away from the area of concern 
(the center of the area), only his lateral position error 
(which is univariate) was considered. 

lt was assumed for this analysis that the lateral 
component of pattern error is nearly normal, and that it can 


somehow be calculated or estimated. 


The expression for the total effort expended in flying 


the spokes of a sector search is 


: of 


pat 1 2 
ra EZ IeineS 
g | | ore 
E= 2Wno fa € rdrdo (31) 
| P 
0 O 


(rsinn = Loterol ronge } 
goxt - Oo Lateral Doecteta error ) 
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The derivative of this expression with respect to r was 
used as the search effort density function for a_ sector 


search with pattern error. 










iy 7 
oy @Ava | ¢ 207° 
j b ane cae 
ex “= OU 
For the analysis” this basic density function was 


computed with sweep width, standard deviation of pattern 


error, and number of spokes all equal to 1. To find the 


density function for any particular set of parameters, the 
density is multiplied by W and n, and the horizontal scale 
is multiplied by the standard deviation of the pattern 


error. 


C. OPTIMUM SECTOR SEARCH 

The plot of this density function is found in Figure 
12. An inspection of this plot will show that it is close to 
a quadratic function out to r = 3 standard deviations. 
Since Koopman's ideal search is a quadratic, some guidance 
was provided as to the optimum radius for the sector search, 
based on the pattern error. It was then necessary to adjust 
the number of spokes to arrive at the _ ideal effort 
expenditure to match this radius. Koopman's ideal search 
effort density is a function of search radius and target's 


position error: 
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Figure 12. Search Effort Density in a Sector Search 


with Pattern Error. 
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Zi Ps 
(tL =f a (Ee 


Clxy) Ad eee 
ae a oe 
and the radius is a_ function of target position error 


and available effort: 
v2 EE] GF 


ae 


Q = 


Solving for effort, E: 


eae 
BS See aE 
Gi 
However, in the case of the sector search, L, the 


number of miles flown, Its simply the number of spokes, n, 
times the radius. The chords flown at the perimeter are 


disregarded in this analysis: 
L=axXon 
Substituting, and solving for n: 
a 
gee [iil 


Ace ay 


Ideally, then, stnce the optimum a Is known to be 
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close to three standard deviations: 


(33) 





or, when all errors are expressed as probable errors, 


Cp = eH Gr 
= 3 
Py hoe Se (34) 
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Tables l “is a 


conditions of 


pattern error. 


total 


tabulation of this function for various 


probable error of position c, and 


The n derived from the table must be divided 


by sweep width W to obtain number of spokes. 


From the 
and W = 
miles, then search 


Interpolating 


Number of spokes = 


and the angle between adjacent search legs = 360/10 = 


degrees. 
Total 


10 X 58.5 


plus the amount of flying at 


fly the pattern, 


somewhat less than 


585 (t5b 


This table for 


accurate means of 


function glven tn 
error, and thus 
use the table. 
if 


Even it 


parameters of a 


example of the sunken 


10.0 miles. 


in the table, 


accuracy for this parameter 


is 


sector 


fishing boat, c = 31.9 

lf the searcher's pattern error were 13 
radius R would equal 58.5 miles. 
n = 100, then 
100/10 = 10 
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search miles flown would be 


585 
the perimeter (in order to 


the aircraft must traverse a distance 


one-half the circumference of the area): 


XT = 585 + 184 = 770 miles. 


sector search will not be useful until an 


predicting pattern error ts found. The 


(3¢ ) is extremely sensitive to pattern 


is required to 


not actually used to determine the 


search, the table does serve to 


68 


SECTOR SEARCH 


TABLE OF OPTIMUM RANTS &ND NUMBER OF SEARCH LEGS 


ERROR, C 
39 35 40 45 50 


PNSITYON 
20 25 
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Divide quantity in "N X W' table by Sweep 
h legs. 


Width to obtain number of searc 


Optimum Sector Search. 


Table 1. 


illustrate that accurate navigation does not necessarily 
produce the optimum sector search, and that only under the 
unusual conditions of an extremely accurate estimate of 
target position is it desireable to reduce pattern error. 


On the contrary, [In most instances, effectiveness may be 


enhanced by increasing pattern error. 
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VI. ONCLUSTONS AND RECOMMENDA NS 


A. THE SEARCH PLANNING PROGRAM 
lt is recommended that the Coast Guard begin immediately 
to make full use of the Search Planning Program now 
operational and undergoing testing at the Fleet Numerical 
Weather Central, Monterey. A guide for users of this program 
is included as Appendix A. 
1. Precautions 
The following cautions are appropriate for users of 
the information provided by this program. 
a. Region of Validity. 

Because the weather data points used in this 
program are approximately 200 miles apart, drift 
calculations provided for areas within 100 miles of shore 
should be double-checked by hand before using. If hand 
calculations produce a larger value for total drift, then 
the hand-calculated answers should be used in lieu of the 
computer solution. 

b. Time of Datum 

Weather predictions in the format reqtined for 

this program are not produced by the Weather Facility for 


more than 4&8 hours in advance. Accordingly, users should 
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avoid requesting plots for time periods beyond 36 hours from 
the time the request is submitted. Ifa plot is requested 
which is later in time than the last computed weather, the 
information will be provided for that plot and a caution 
statement will appear in the message reply. 
c. Numbers of Initial and Datum Positions 
No more than nine initial positions may be 
entered, and no more- than nine datum times will be 
calculated, in any one run of the program. 
The user iS encouraged to request whatever 
number of datum points may be useful to him. Thus if a 
search is planned for 0800Z, the search planner may wish to 
request additional data for target position at 1000Z, 1200Z, 
1400Z, and 1600Z. 
2. Leeway 
it is suggested that further experimentation be done 
in the area of survivor craft leeway to determine’ the 
general shape of the leeway function. It is expected that 
this shape will be close to that of the liferaft drift 


tables, and to the formula of Witting given in Ref.3. 


k \| Wind (35) 


and which is used in the time-sharing version of the 


— 
i 


program for all leeway. This same function is used in both 


versions to approximate liferaft leeway. 
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When better data is available, the program at Fleet 
Numerical can be changed to accomodate the curved leeway 
function. 

3. JLesting 

It is suggested that a rigorous program of testing 
and evaluation of the results of the Search Planning Program 
be initiated. If the program proves to be more accurate 
than other methods, its contribution to the search planning 
process, and to oceanographic research in general, will be 
substantial. 

4. Qperational Vse in the Coast Guard 
Although the program will be available to the Coast 
Guard via the facilities of the Fleet Numerical Weather 
Central, there are further possible advantages of such 
program which are not available through these facilities. 
a. Optimum Search Areas 
The computer method of solving for optimum 
searches shown in. part IV of this thesis could be added to 
the Search Planning Program. 
b. Other Applications 
Computation of on-scene aircraft endurance, 
time enroute search areas, division of search area among 
different search units, maximum rescue coverage intercepts, 
fixes from radio direction finder bearings, and many other 
calculations familiar to search planners could all. be 


available instantly from a remote computer terminal. 
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The Coast Guard Amver computer center at Governor's 
island New York regularly receives weather data by wire from 
Fleet Numerical Weather Central which is used to produce 


weather maps. The fields required for the Search Planning 


Program could also be transmitted to a central Coast Guard 
computer facility. Remote terminals connected to this 
facility and tInstalled in Rescue Coordination Centers 
(RCC's) could provide the output from this program, Amver 
surface picture plots, and other useful data directly to 


search planners. 


8, DETERMINATION OF SEARCH AREAS 

it is recommended that the findings of parts IV and V of 
this thesis be evaluated by the Coast Guard and applicable 
portions be included in the National Search and Rescue 


Manual. 


C. FUTURE RESEARCH 
There are many unexplored areas [tn the field of Search 
Planning in which further research may yield valuable 
results. 
ren ee in 
Although the type of point-datum search presented in 
parts IV and V_ are computed in two dimensions, there are 
times when the target's distribution is centered ona line 


of position (or trackline), which presents a- problem 
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essentially one dimensional in nature. 


pe 


if a boat departs point A intending to follow a 
certain track AB to point B, and does not arrive at point B, 
and if no other information is available as to his probable 
position, a one-dimensional or trackline search problem is 
presented. 
a. Partial Solution 
Assuming that it is decided to begin the 
searching efforts by executing just one sweep along. the 
length of the trackline, what is the probability that the 
target will be located? The following analysis is for the 
inverse cube law of detection. 
Given the standard deviation of the target's 


error of position ©, and given sweep width, W, where 


ee 2a ez 
(36) 


Solving for m, which includes visibility and 


other factors, 


~w) 





Ww 
Aa an Q TT 
(37) 
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The probability of detecting a uniformly 





distributed target is aK 
oe 
OC. 
We (38) 
fe 2 
a nagar 417 x 
(39) 
Introducing a new variable, 6=-= 
oo 
a 2 i 
Leer ie 
(40) 


The probability that the target is located at 


some distance x from the trackline may be assumed to be 


{ ees 
x)= —— Ney 
P ) eae Ooms © (41) 


Thus, the probability of detecting this target 
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which has the solution 
_ ¥2- W 


— (45) 


b. General Solution 
The solution for probability P of detecting a 
target on a trackline search for more than one sweep 
Involves a more complicated expression, one which does not 
lend itself to such a_ simple solution. Since a trackline 
search Is executed one sweep at a time, a general solution 
for the optimum distance of each sweep from the trackline is 
Indicated, and should be the subject of further study. 
2. Qptimum Search Speed 
The Search and Rescue Manual prescribes a "moderate" 


speed for searching. Koopman's formula for sweep width, from 


(9) gives 
ZkKh 
eS Us 
This Indicates that sweep width Is Inversely 


aaport onal to the square root of the relative speed of the 
target, w. When search speed Is very large in comparison to 
target speed, w takes on the value of the searcher's speed. 
As w Increases, however, the amount of search path length L 
that can be covered in a_ certain period of time increases 
proportionally: 


| = K, Ww 
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It follows that the amount of search effort, Wx L, 
that can be effected by a certain number of aircraft ina 


fixed period of time Increases as the square root of search 








speed: 
W_ . x = INL 
S ‘ A 
— WE 
f\ = Ko 
\NL Az yah ew ees (47) 
Kz ia K2 


lt may be argued, then, that when time ts a factor, 
such as where persons In distress may not’ survive for long 
periods of time, searching should be accomplished at’ the 
highest speed practicable. In addition, since the coverage 
of a designated search area will take less time at a higher 
speed, lookouts may encounter less’ fatigue, and the 


effectiveness of the search may be enhanced. 


D. CONCLUSION 
In concluston, there Is a need for more. scientific 
research In the practical aspects of search planning, and a 


corresponding need for the education of search planners in 


the fields of probability and search theory. 
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APPENDIX A 


USER'S GUIDE TO 
SEARCH PLANNING PROGRAM 


1. Description 

The MONTEREY SEARCH PLANNING PROGRAM {fs a_ computer 
program which solves open ocean drift plotting problems, and 
provides datum, minimax datum, Total Probable Error of 
Position, and Search Radius for both simple and complex 
searches. The use of this program is available upon 
request, on atest basis, to all Rescue Coordination Centers 
through the facilities of the U.S. Navy Fleet Numerical 


Weather Facility, Monterey, California. 


a. Leeway Code 
For this program, leeway is expressed as a two-digit 


code as shown in Table l. 
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00 
10 


leeway Code 


No Leeway 

Liferaft with Drogue 

Liferaft without Drogue 
Liferaft (Minimax) 

Downwind drift (minimax) 
Minimum = 1% of wind velocity 
Maximum = Second digit as % of 


Wind Velocity. 


20,30,etc. Downwind drift at 2%,3%, etc. 


/ i Oa Ie 


of the wind velocity 

Minimax drift at some fixed percent 
of wind velocity equal to first digit; 
Deflection to the right and left 

of the wind line equal to 


10 degrees times the second digit. 


Table 1. Leeway Code 


Example: G3 = 4% of wind velocity 
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+ 30 degrees off wind line 


b. Number of Starting Points 
The number of starting (or distress) positions must 
be indicated. The maximum number that can be supplied is 
nine. 
c. Number of Datum Times 
The number of datum (or search) times for which 
plots are requested must be supplied. The maximum number of 
these is also nine. 
d. Starting Point Data 
Starting positions must be entered in chronological] 
order and each point must be fully identified by: 
(1) Date-time group (time zone Zulu) 
(2) Month 
(3) Latitude and Longitude in degrees and minutes. 
Tenths of minutes must be specified. The number’ of these 
starting positions must agree with paragraph 2.a. 
e. Datum Times 
Complete Date-time group Cincluding month) for each 
datum time must be supplied, in chronological order. The 
number of these datum times must agree with paragraph 2.c. 
f. Last Known Position 
This is the last known position of the search object 


(see Para. 611, National Search and Rescue Manual) to be 
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used in computing DR position error. This position must be 
complete with degrees, minutes, and tenths of minutes. If 
DR position error ts not applicable, enter one of the 
initial positions from 2.d. 
g. Error Factors 
These factors must all be specified to at least one 
decimal place, in the following order: 
(1) Probable navigation error of search craft. 
2) Initial probable error in position of the 
search object. | 
(3) Percent of DR distance to be applied as ODOR 
error (See Para. 611, National Search and Rescue Manual). 
Enter 0.0 if not applicable. 
h. Search Designators 
These should be one-digit numbers, from 0 to 5, 
which describe which Safety Factors are to be applied to the 
Total Probable Error of position for each datum time. One 


digit must be supplied for each datum time. 


0 No search radius is to be computed for 
that datum time 
1, 2,357495 The Safety Factor corresponding to 


search 1,2,3,4, or 5 is to be applied. 


3. Message Request Format 


A message request to Fleet Numerical Weather Central 
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will insure that the data is entered correctly into the 


computer: 


FM CEA 
TO FLENUMWEACEN MONTEREY 
BT 
UNCLAS 
DISTRESS HOMEY 65 OVERDUE 
1. REQUEST SAR PLOT. 
2. DATA: 
Ay 12 eee 
B. §@39989Z AUG,3529.9N,85216.5W/93129HZ AUG, 
3680.9N, 95148. 9W// 
C., 0316882 AUG/$322H98Z AUG/B4P2HBZ2 AUG// 
D. 3599.8N,85149.9W/5.8/18.8/8.18/192// 


This sample message requests a minimax plot of a 
liferaft from two initial positions, with data requested for 
three different times (para. 2.A.). The Initial positions 
and date-time groups are given in para. 2.B. (note that 2 
digits are used to express degrees latitude and 3 digits are 
used for degrees longitude). 

Paragraph 2.C. gives the date-time groups of the datum 
times. Paragraph 2.D. gives the last known position, error 
factors, and search designators. for this example, the 


first search is to be computed for 031600Z anda ee second 
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search for 040200Z. No search radius is to be computed for 


0322002; positions and position errors will still. be 


provided. 


GW, Caution: Offshore Regions 

This program is designed primarily for open ocean 
search; when data is supplied for water areas within 100 
miles from shore, answers’~ should be checked by hand 
calculations. Requests for SAR Plots within 10 miles from 
shore should be avoided. 

Weather information used in this program Is applicable 


to the Northern Hemisphere only. 
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